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We show that a string-inspired Planck scale modification of general relativity can have observable 
cosmological effects. Specifically, we present a complete analysis of the inflationary perturbation 
spectrum produced by a phenomenological Lagrangian that has a standard form on large scales but 
incorporates a string-inspired short distance cutoff, and find a deviation from the standard result. We 
use the de Sitter calculation as the basis of a qualitative analysis of other inflationary backgrounds, 
arguing that in these cases the cutoff could have a more pronounced effect, changing the shape 
of the spectrum. Moreover, the computational approach developed here can be used to provide 
unambiguous calculations of the perturbation spectrum in other heuristic models that modify trans- 
Planckian physics and thereby determine their impact on the inflationary perturbation spectrum. 
Finally, we argue that this model may provide an exception to constraints, recently proposed by 
Tanaka and Starobinsky, on the ability of Planck-scale physics to modify the cosmological spectrum. 



I. INTRODUCTION 

In the last two decades, particle physics has become 
an indispensable part of cosmology. In fact, one of the 
strong motivations for studying particle physics theories 
that go beyond the standard model and incorporate grav- 
ity is that they may shed light on the nature of the cos- 
mological singularities arising in general relativity. There 
is "widespread hope that in one form or another these and 
other cosmological considerations may one day allow us 
to test physical theories whose fundamental scales are 
now, and perhaps forever, beyond the reach of conven- 
tional accelerator experiments. 

Inflationary cosmological models, in particular, signif- 
icantly highlight the important role of microphysics. For 
example, since the pioneering work of Q it has been 
known that quantum field fluctuations in the early uni- 
verse are stretched by inflationary expansion to scales of 
astrophysical relevance, providing a gratifying first prin- 
ciples mechanism for structure formation. Galaxies, from 
this viewpoint, are quantum fluctuations writ large. 

In order to solve the standard cosmological puzzles, a 
minimum of 60 e-folds of inflationary expansion must be 
invoked, but in many models this number can be much 
larger. Taking this at face value, it means that today's 
scales of cosmological relevance expanded from Planckian 
or sub-Planckian scales at the onset of inflation. Inflation 
may therefore provide a kind of Planck scale microscope, 
stretching the smallest of distance scales to observably 
large size. 

It is possible, however, that in the process of such in- 
flationary expansion the effect of Planck scale physics 
gets washed out, being diluted by the very growth of 
scales which potentially makes them visible. A similar 



phenomenon has been observed in black hole physics. 
If one traces the history of a Hawking radiated photon, 
one finds that it gets ever more blue-shifted toward the 
moment of its emission, and hence one might think it 
could carry an imprint of extremely high energy physi- 
cal processes. In reality though, a number of studies 
have concluded that such short distance physics does not 
have any impact on the low energy features of Hawking 
radiation; heuristically, short distance modifications are 
washed out by the memoryless process of thcrmalization. 

In a cosmological context, the situation in this regard 
has been less clear. Brandenberger and Martin ||] car- 
ried out a study of the impact of various models of trans- 
Planckian physics (in fact, the same models considered in 
the black hole studies just mentioned) on the spectrum of 
density perturbations in power law inflationary models. 
They found that while some hypothesized Planck scale 
modifications to ordinary field theory yield no late time 
consequences (similar to the black hole conclusion) some 
do, indicating a cosmological sensitivity to short distance 
physics. On the other hand, it was argued in [|| that 
there should be no change in the perturbation spectrum 
if the proposed modifications still yield the adiabatic vac- 
uumQ. 

In this paper we take up the issue of cosmological sen- 
sitivity to short distance physics, but from a different 
approach. Namely, rather than considering the ad hoc 
short scale modifications studied in , we focus on an 



*Note that the definitions of adiabaticity used in M and 
j^] are slightly different, leaving the situation somewhat 
ambiguous. 
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interesting model introduced in M , whose short distance 
physics is designed to incorporate a minimum length, 
a development naturally inspired by string theory. We 
present a full analysis of the perturbation spectrum pro- 
duced by the presence of such a minimum length in a 
de Sitter background. As with any calculation of the 
inflationary perturbation spectrum, the modes are nor- 
malized well inside the horizon, and the calculational task 
is to track their evolution until they are far outside the 
horizon, as the perturbation spectrum is fully determined 
by their asymptotic values. We solve the mode equation 
numerically to avoid invoking any new analytical approx- 
imations, and fix the initial conditions for the numerical 
solution by matching to an approximate analytical so- 
lution of the mode equations that is valid at very early 
times. 

While completing this work, Kempf and Niemeyer |?|] 
posted a complementary calculation of the spectrum in 
the presence of a short distance cutoff. We contrast our 
approach with theirs, showing that the numerical evalu- 
ation of the mode functions leads to a complete descrip- 
tion of the spectrum. We also expect that the approach 
used here will generalize to the broader problem of de- 
termining the impact of trans-Planckian physics on the 
inflationary perturbation spectrum, which can manifest 
itself as a change to the dispersion relation [||, or the 
initial normalization of the modes or choice of vacuum. 

In what follows, we briefly review the construction of 
a model with a minimum length, and discuss the ana- 
lytical properties of the perturbation equation. In par- 
ticular, following H we show that since there is a min- 
imum length, each mode is "created" at a finite time in 
the past, defined by the moment when its physical wave- 
length first exceeds the minimum length. We then show 
how to extract an approximate analytical solution that 
holds near this early time. This allows us to "match" 
this solution to a numerical solution of the full mode 
equation, and compute the amplitude of the perturba- 
tion spectrum as a function of the minimum length. We 
find that, indeed, there is an imprint of the short dis- 
tance cutoff on the perturbation spectrum. We also con- 
firm that the scenario we study does not violate bounds 
coming from late time particle production and quantum 
mechanical self consistency. More generally, we empha- 
size that the effect of any short distance modification can 
be encoded in two wavelength dependent functions, pro- 
viding a phenomenological approach for systematically 
analysing deviations from the standard prediction of the 
power spectrum. We conclude with a brief discussion of 
the spectrum we might expect in a more general space- 
time, and a summary of our conclusions. 

We stress at the outset that our intent is to show ex- 
plicitly that a modification to conventional physics at the 
Planck scale can have an observable effect on precision 
cosmological measurements. While of interest in its own 
right, the model M on which we focus should be viewed as 



one concrete example in which such calculations can be 
reliably performed, allowing us to establish definitively 
that cosmological observations may be a window onto 
Planck scale physics. 

II. PERTURBATIONS WITH A MINIMUM 
LENGTH 

Motivated by the stringy uncertainty principle j5|, a 
phenomenological Lagrangian which incorporates a short 
distance cutoff was recently proposed by Kempf J(|. In 
this approach, the short distance cutoff is modeled by 
modifying the usual commutation relation to 

[x, P ] =ih(i + p P 2 ). (l) 

The parameter (3 is related to the minimum distance 
Ax min by Ax rnin — y/]3. A Lagrangian suggested by 
Eq.(|l]) was discussed in ||. In this model, the tensor 
mode Uj, obeys the following equation of motion: 

„ v' ( a" a' v'\ 
K v K \ a a v J 

where a denotes the scale factor and the prime denotes 
differentiation with respect to conformal time r\. Our u^. 
is equal to a 2 4>^ from while k % = ap l e~^ p2 1 2 where p l 
is the Fourier transform of the physical coordinates x l , 
and 

a 2 p 2 e^P 2 

The cutoff is defined by requiring that p 2 < 1/(3, and 
is motivated by the minimum distance in string theory. 
Tensor perturbations with different comoving wavenum- 
ber k reach the cutoff p 2 = 1/(3 at different confor- 
mal time ?/fc, where o?{rjk) — {3k 2 . For de Sitter space, 

r? fe = -l/y/e/3H 2 k 2 . 

Note that in evaluating the derivatives of Uf. with re- 
spect to 77, we are holding k (and not the usual comov- 
ing momentum k) fixed with time. It is therefore conve- 
nient to express p and v in terms of k by introducing the 
product-log, or Lambert W function Q , which is defined 
so that W(xe x ) — x: 

a 2 W{z) z/ _ a/_W(z) (5 + 3W(z)) 

fi ~~J (i + w(z)f ' (\ + w{z)Y ■ 

(4) 

where z — —(3k /a 2 . The W function has an essential 
singularity when its argument is equal to — 1/e, and this 
corresponds to the precise moment (for a given k) when 
f] = rjk- Let us examine in more detail the nature of the 
singularity in the equation of motion at 77 = r/k- The 
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Lambert W function W(z) has a series expansion near 
the branch point z — — 1/e ||: 



W(z) 
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(5) 



where p = \j2iez + 1). The series converges for |p| < y2- 
The singular point atr) = rjk is irregular because the coef- 
ficients of 7t~ and are not analytic in r\ — rjk- However, 
the non- analytic piece is less singular than 1/(77 — 77/.) and 
is therefore subdominant. We can solve for the leading 
behavior of uz by extracting the most singular terms of 
the equation of motion. First, write rj = 77/5(1 — y), so 
that the leading terms in the equation of motion give 
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where dot denotes derivatives with respect to y. Written 
in this form, the ii^ and ut terms both have divergences 
that scale as 1/(77 — 77^), so the equation for the leading 
behavior of ut actually has a regular singular point at r/k- 
Since any second order differential equation has a 
power-law expansion about a regular singular point |[o[ , 
we can explicitly construct the two independent solutions 
of Eq. (^) in the vicinity of r/k ■ Let take the form: 



, c n y ' 



(7) 



The indicial equation is obtained by inserting this ex- 
pression into into Eq.lQ), expanding about y — and 
indentifying the lowest order terms, which are propor- 
tional to y a ~ 2 ■ This gives a — or 3/2. Therefore, for 
small y, 



ci + c 2 y 



3/2 



(8) 



The coefficients are constrained by the Wronskian 
condition which follows from the canonical commutation 
relation \(f>z, TTf] = i5 3 (k — r) for <j)f. — uy 'a and its con- 
jugate momentum TTf HQ] [J} 



ujfaKfa) - u t(vH(v) = 1 (1 - PP 2 ) exp(--/3p 2 ). 



(9) 



Hence, 



tin the model of the commutation relations of x and p 
are modified as in Eq. (|l|) , but the field commutation relations 
are not changed from their standard form. While one can seek 
a particular interpretation, here we simply view this as pro- 
viding an interesting, short distance modification of the field 
equations that allows for a reliable calculation of cosmological 
implications. 
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c x ci = iri k -e 1 . 



(10) 



The Wronskian condition mixes the two power series so- 
lutions with complex coefficients. A priori, there is a one 
parameter family of solutions of C1.2, corresponding to 
different choices of the vacuum state. In the usual infla- 
tionary scenario, the natural choice of the vacuum state 
is the so-called Bunch-Davies vacuum which reduces to 
the Minkowskian vacuum for wavelengths much shorter 
than the Hubble scale. We now motivate a natural choice 
of vacuum in the present context by comparing with the 
Bunch-Davies vacuum. Recall that in the standard set- 
ting, the high frequency limit of the mode equation takes 
the form 



K(v) + ^KvWiv) = 0, 



(ii) 



whose solution is well approximated by the WKB form 



1 r 



(12) 



<C 1 are 



if the adiabatic conditions — *- -C 1 and 
both satisfied. The Bunch-Davies vacuum amounts to 
choosing the Heisenberg form of the field operator 71(77, x ) 
to be 



t2(?7,x) 



d 3 k 
(2tt) 3 / 2 



(u^a k e 



ikx 



(13) 



with a vacuum state |0) satisfying afe|0) = 0. In our case, 
if we could ignore the in (^) , similar reasoning suggests 



where, from 
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exp(-i / u> k (r)')dri'), 



we see that for 77 
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and 
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(14) 



(15) 



(16) 



Af3H 2 2' 

This would give a Bunch-Davies-like vacuum of the form 



\AA 



2 \ !/ 4 



exp(— 2iy/ Ay) 



(17) 



Of course, though, we can not ignore the uj. in (^|), and 
hence this does not yield a solution to our modified mode 
equation. However, by modifying the prefactor of the 
exponential, we can in fact construct a solution of this 
form, 



f(v) = 



+ iA^/y exp(-2»yj4y). 



(18) 
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This solves Eq. (^|); equivalently, in the series expansion 
of F(y) the y 1 ! 2 term cancels out, and F(y) is seen to 
be a linear combination of the two power series solutions 
found above with ci = ^ and c 2 = — i|A 2 . 

The general solution is a linear combination of the pos- 
itive and the negative frequency modes 



u- k {y) = C + F{y) + C^F*{y) 



(19) 



with the constants C± constrained by the Wronskian con- 
dition 



IC+I 2 
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y^fctf IP 2 2 
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(20) 



Comparison of the osciallatory part of F(y) with that of 
the Bunch-Davies vacuum suggests that C_ = 0. How- 
ever, the normalization of F(y) is not 1/v2cj since the 
adiabatic condition is not satisfied. To be specific, 



,2 



(21) 



Therefore for y ~ 0, no matter how large ^4 is (or, more 
importantly, how small (3 is), the adiabatic condition 

|^%| << 1 is violated. Thus these initial conditions are 
not smoothly connected to the "standard" form of the 
mode equation when it is well inside the horizon, even 
in the limit (3 — > although the Lagrangian is smooth in 
the same limit. The discrepancy arises because whenever 
(3^0, the mode's evolution begins at a finite confor- 
mal time. Moreover, we expanded the Lambert function 
W(z) with z — —(3k 2 /a 2 around z — — 1/e, and this ex- 
pansion is only convergent when z < 0. If (3 = 0, the 
argument of the W function is zero, and this expansion 
cannot be used. 

We have found the leading behavior of ut around rj ~ 
r/fe. The equation of motion for the tensor mode is 
solved only up to order 1/y. The residual terms can still 
be significant for rj ~ rjk- We deduce the subleading 
bahavior of by the method of dominant balance [ pd| . 
Define uj,{y) = F(y) (1 + ei(t/)) and extract the most 
singular terms in the equation of motion for e± : 



e i - ;r ei 



3A- 3/2 



o. 



(22) 



which gives 



e l {y) = (2A-l)y 3 / 2 [logy- 



This is indeed small compared with the leading term for 
small y. With this correction, the equation is solved up 
to order 1/y/y, but there are still residual ln(y) terms. 
The solution is further improved by the next subleading 
order e 2 (y), where = F(y)(l + ei(y))(l + £2(2/))- The 
solution for 62(2/) is 
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FIG. 1. The top figure shows the spectrum (P 1 / 2 ) as a 
function of /3 with C-/C+ = —.5, while the lower plot shows 
the spectrum when C_ = 0. It is only in the latter case 
the computed value of the spectrum smoothly approaches the 
usual value in de Sitter space, In this plot H = .1, and with 
(5 = we would expect P 1/2 = 0.0159155. 



€ 2 = — (105 - 33CL4 - 112iA 3/2 j y 2 
7 -(2A-l)y 2 \ogy. 



(24) 



The equation of motion for is therefore solved up to 
terms that vanish as 7/ — > 0. 

When the mode is well outside the horizon, p <C H, 
and see that ii(rj, p) — ► k 2 and 1/(77, P) ^ 1- Consequently, 
Uk(il) ~ a (v)t which reproduces the standard late time 
limit in the usual case with (3 = 0. This limit defines the 
power spectrum, 



(25) 



(23) which is evaluated when u/- is well outside the horizon. 



III. NUMERICAL RESULTS 

We match the analytical form of the solution valid 
when the mode is well inside the horizon to a numeri- 
cal evaluation of the full mode equation. In principle we 
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FIG. 2. The p dependence of the spectrum, P k ' , is plotted 
against f3, with H = 0.1. For a de Sitter background, the 
spectrum is independent of k. 

might hope to evaluate the mode equation numerically 
from an initial time of y = 0, but the coefficients of uj, 
and u'r are infinite at this point so it is easier to match 
the numerical solutions to the approximate analytical so- 
lutions at yo, a small but finite value of y. By varying 
yo we can ensure that our results do not depend on our 
choice of starting point, and that we are therefore "close 
enough" to y = 0. 

A subtle point that arises during the integration is that 
the subdominant terms in uz contribute terms of order 
y 3 / 2 log(y), which contribute y/ylogy in the derivative, 
it ~ . While these terms do go to zero in the limit where yo 
is vanishingly small, they approach zero comparatively 
slowly. Since we are using a small but not infinitesi- 
mal value of yo j we match the numerical solution to the 
"corrected" asymptotic solution, which includes the low- 
est order logarithmic terms. We evolve the mode equa- 
tions numerically using the Bulrisch-Stoer integrator im- 
plemented in Fortran, and from the asymptotic values of 
|tifc/a| we obtain the spectrum, P 1 / 2 (fc). 

In general the spectrum is ^-dependent, but for the 
special case of de Sitter inflation, the background space- 
time is time translation invariant, which implies that the 
spectrum should not depend on k (which is equivalent to 
k at late times). This is manifest from the mode equation 
and, in practice, we verified the code by solving the mode 
equation for multiple values of k and found that the nu- 
merically computed spectrum was indeed scale invariant 
to better than 1 part in 10 6 . 

In the previous section, we saw that imposing the 
Wronskian constraint led to a one parameter family of so- 
lutions, but that a comparison with the standard Bunch- 
Davies vacuum suggests the choice C_ = 0. We begin 
by analysing the consequences of relaxing this choice. In 
Fig. (Q) we plot the spectrum computed for small val- 
ues of (3 with C_ = 0, and compare this plot to a cal- 
culation with a finite value of C_. In the former case, 
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FIG. 3. The values of D+ (top) and D- are extracted for 
the spectrum shown in Fig. 1 with C- — 0. As f3 — > 0, 
D+ — > 1 and D- — > 0, confirming that as the Hubble length 
becomes much larger than the minimum length the spectrum 
approaches the standard result. 

the spectrum smoothly approaches the de Sitter result 
P 1 / 2 — H/2n, but does not do so for any finite value of 
CL. 

Fig. |2| displays the (3 dependence of the spectrum. For 
moderate values of /?, the spectrum can be either lower 
or higher than the "standard" value, but for large /?, 
P 1 / 2 approaches zero. For sufficiently large (3, the overall 
amplitude of the perturbation spectrum is significantly 
reduced. In de Sitter space the power spectrum is k- 
independent, so the only effect is an overall normalization 
dependence on (3 through the dimensionless combination 
\/]3H. Assuming one has an independent measure of H, 
we see that the short distance cutoff specified by (3 does 
indeed leave an imprint on the spectrum. In more real- 
istic inflationary models, as we indicate later, we believe 
the imprint will be /c-dependent and hence also affect the 
shape of the power spectrum. 

IV. MODE MATCHING AFTER HORIZON 
CROSSING 

The physical origin of the change in normalization of 
the power spectrum can be simply understood in terms 
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of the asymptotic behavior of the mode solutions. In the 
long wavelength limit, the mode equation asymptotes to 
the standard mode equation for cosmological perturba- 
tions. Therefore, regardless of the nature of the short- 
distance physics, the mode function takes the general 
form: 



1 



u k (v) = [D+H v (- 



kr,) + D_H*(-kr,)], (26) 



where H v is a Hankel function of the first kind, and v = 
3/2 in de Sitter space. The constants D± are constrained 
by the standard Wronskian condition arising from the 
canonical commutation relations for the field operator: 

„ duk dut 



drj k dr] 



(27) 



With this constraint, the form of the mode function is 
completely determined (up to a phase) by the selection 
of a vacuum. The standard case, that of a Bunch-Davies 
vacuum at short wavelength, corresponds to the selection 



£>_ = and D, 



1, with our choice of normalization. 



\D+ + D. 



The key attribute of this construction is that regard- 
less of the nature of the short-distance physics, the long- 
wavelength behavior of the modes is completely encoded 
in the coefficients D± . In general the D± are fc-dependent 
but in de Sitter space the situation is particularly simple, 
since the D± are constants. 

In the standard analysis, the coefficients C± that deter- 
mine the vacuum at short wavelengths are identical to the 
D± found by matching the mode functions to the long- 
wavelength limit, Eq. (p6|). Modifying the short distance 
physics can break this correspondence, and the choice of 
vacuum C_ = in the short wavelength limit can corre- 
spond to a phase rotation in the long wavelength limit, 
with D_ =/= 0. This changes the normalization of Mfc, since 
the power spectrum depends on the integration constants 
as 

fc 3 

= 2^ 

in the long wavelength (fc — > 0) limit. In de Sitter space, 
the D± are constant, so the rotation alters the normal- 
ization of the power spectrum, but does not introduce 
any k dependence. However, if the background deviates 
from de Sitter the D± are, in general, fc-dependent, which 
will alter both the shape and normalization of the power 
spectrum as well as the short distance physics - perhaps 
in a dramatic way Naturally, if a modification to 
short-distance physics alters D±, the perturbation spec- 
trum will change. This will lead to constraints on the 
form and magnitude of the modifications to D±. 



V. PARTICLE PRODUCTION BOUNDS 

In addition to the perturbation spectrum, other ar- 
guments place cosmological constraints on modifications 



(28) 



to short-distance physics. For example, Starobinsky H] 
has recently shown on quite general grounds that even a 
very tiny rotation away from the Bunch-Davies vacuum 
in the current universe would result in unacceptable pro- 
duction of relativistic particles. The Starobinsky bound 
corresponds to 



\D. 



< 



(29) 



where the coefficient D_ is evaluated in the limit of a 
scale much larger than the minimum distance y/]3 but 
still much smaller than the current horizon size Hq 1 . 
However, this bound applies to the cosmic vacuum at 
late times, not during inflation. If \D-\ is dependent on 
H, it can be significant during inflation, thus affecting 
the power spectrum, but small enough universe today to 
satisfy the Starobinsky bound. 

As it happens, an H dependent value of £>_ is precisely 
what the model considered in this paper predicts. The 
magnitude of the rotation depends not on the absolute 
physical scale of the cuttoff, but on the ratio of the cutoff 
scale to the horizon size, \fj3~H. If the horizon size and 
the minimum length are comparable during inflation, the 
power spectrum is significantly affected, but the rate of 
particle production today will be strongly suppressed, by 
the much lower value of H in the present universe. Fig. 3 
shows the dependence of the coefficient D_ on the ratio 
f3H 2 in the numerical solution of the mode equation. The 
numerical results imply \D- \ oc j3 n , with 0.45 < n < 0.5 
when (3H 2 <C 1. In the "real" inflationary universe, we 
hold (3 fixed and reduce H, so this result implies that as 
the universe expands the D + term will become negligible. 

On the basis of our present calculations, the exact de- 
pendence of D- on j3 remains unclear, since the value of 
n appears to be very weakly dependent on (3: if we eval- 
uate n over a few decades in {3 (for fixed H) n appears to 
slowly approach 0.5 as this range is moved to smaller and 
smaller j3. Thus, it is not unreasonable to assume that 
\D—\ oc /3' 5 in the present universe, and that the particle 
production rate today is consistent with the Starobinsky 
bound. 

We do not pursue this in more detail as the current 
model is simply a heuristic construction that mimics 
what we might expect from a more rigorous string theo- 
retic description of spacetime. However, we can conclude 
that it is plausible that a minimum length leads to a 
time dependent modification of the backgroun in a uni- 
verse where H is not constant, and that the value of D_ 
will be time dependent making the particule production 
bound much less onerous. 

Also of relevance is Tanaka's analysis pTq] , showing 
that particle production during inflation can significantly 
change the background evolution of the spacetime. The 
perturbation to the stress-energy during inflation can be 
expressed as: 
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k\D- 



_fk_ 

d 3 P , n ,2 

TP \D-\ 

(2^) 3 ' 



\D-\ 

P 2 



(30) 



where p is the physical momentum p = (k/a), with the 
momentum cutoff p < (l/^/JT). If perturbation theory is 
to be consistent, the contribution from the stress-energy 
from particle production must be subdominant, 
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Tanaka's analysis assumed \D_ \ 
on f3 of 

1 



(3 



« y/M Pl H 2 . 



(31) 

1 , resulting in a bound 
(32) 



In this model, D— is /3-dependent, weakening of the 
Tanaka bound, and for D_ ~ y/^H it becomes 



pi- 



(33) 



This is easily satisfied if \f]3 is within an order of mag- 
nitude or two of the Planck length. Therefore particle 
production places no significant constraints on the via- 
bility of this model. 



VI. DISCUSSION 

We have presented an accurate calculation of the per- 
turbation spectrum predicted for a de Sitter universe 
where the physics has been modified to include a min- 
imum length. The spectrum is scale independent, which 
follows from the time translation invariance of a de Sit- 
ter universe. This calculation significantly extends the 
contemporaneous qualitative approach of (7) , allowing us 
to determine how the perturbation spectrum is modified. 
Furthermore, it lays the groundwork for accurate calcu- 
lations of the spectrum produced by other models where 
trans-Planckian physics has modified the "standard" evo- 
lution of cosmological perturbations, either by altering 
the dispersion relationship in the evolution equations, or 
by changing the initial conditions. 

If the minimum length is much smaller than the Hub- 
ble length, its introduction has no detectable effect on 
the spectrum. However if these two lengths are within 
an order of magnitude or two, it is possible for the result- 
ing spectrum to differ appreciably from the (3 — limit 
and for various other particle production constraints, dis- 
cussed above, to be satisfied. Various values for the string 
scale and the Hubble scale arise in recent string theoretic 



approaches to cosmology (see for example ]lq-|22]|) and 
hence the effects studied here can potentially be signifi- 
cant in these and other models. 

In de Sitter space, the ratio between the minimum 
length and the physical horizon size is constant. In al- 
most all other inflationary backgrounds, the expansion 
rate is slower than exponential, and the physical horizon 
size will increase relative to the minimum length scale. 
In this case our analysis of the de Sitter background sug- 
gests that the amplitude of the longest modes (produced 
earliest in inflation) will be modified. Shorter modes will 
leave the horizon at a time when the horizon length is 
much larger than the physical cut-off length and their 
amplitude will be unaffected. We plan to return to this 
problem in future work, but our tentative conclusion is 
that the spectrum of primordial perturbations could be 
altered at very long wavelengths by the existence of a 
minimum length scale. Whether this is observable in 
practice will depend crucially on the number of e-folds 
of inflation preceding the creation of the modes which 
are responsible for large scale structure in the observable 
universe. 

We close with a few observations, some rather specu- 
lative, about further studies we plan to undertake based 
on the present work || . 

• All of the calculations in this paper have focused on 
tensor perturbations. It would be worthwhile to extend 
the analysis to scalar perturbations. 

• We find it particularly interesting that since short 
scale modifications yield mode equations that are asymp- 
totic to the standard form on large scales (that is what 
is meant by a short scale modification) the effect on the 
power spectrum of any new short distance physics can 
be encoded in the /c-dependent coefficients D + and D_. 
This provides a phenomenological framework for system- 
atically seeking signals of - and establishing constraints 
on - short scale deviations from conventional physics. 

• For very large j3, the power spectrum becomes arbi- 
trarily small. It is tempting to interpret this as a mech- 
anism for solving the fine tuning problem endemic to 
the inflationary generation of the primordial perturba- 
tion spectrum. A mechanism ensuring that the overall 
normalization of the spectrum is small enough to satisfy 
observational constraints from Large Scale Structure and 
the microwave background will solve this problem. For 
large /?, this model appears to do just that, but we cau- 
tion that this conclusion is likely naive since it requires 
the minimum physical length to be much larger than the 
horizon volume, and we can not be sure the model has 
any physical validity in this regime. Further study along 
these lines may reveal a trustworthy suppression mecha- 
nism. 

• An interesting - but highly speculative - possibility 
is that particle production induced by a finite value of 
D- may both be small enough to satisfy Starobinsky's 
bound, but large enough to lead to new physics. Epochs 
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where this would be particularly intriguing are immedi- 
ately after inflation, where particle production from the 
vacuum is a potential mechanism for reheating the uni- 
verse, and in the present universe if the particle produc- 
tion was efficient enough to modify the equation of state. 
To explain the latter possibility, the equation of state 
determines the relationship between the density and the 
scale factor as the universe expands and particle produc- 
tion reduces the rate at which the density of the uni- 
verse drops with increasing volume. In general, the more 
weakly the density of a perfect fluid depends on volume, 
the more rapidly the universe will expand. Consequently, 
if particle production is efficient enough to alter the ex- 
pansion rate of the universe, it will be increased relative 
to that of a universe with _D_ = 0. This is particularly 
interesting in the light of observational evidence for dark 
energy, which is betrayed by a too-rapid expansion of 
the spacetime background. For an alternative approach 
of understanding the origin of dark energy from trans- 
Planckian physics, see f23||. 
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